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We study a spinless hardcore boson model on checkerboard lattices by Green function Monte Carlo 
method. At half filling, the ground state energy is obtained up to 28 x 28 lattice and extrapolated 
to infinite size, the staggered pseudospin magnetization is found to vanish in the thermodynamic 
limit. Thus the (n,n) charge order is absent in this system. Away from half filling, two defects 
induced by each hole (particle) may carry fractional charge (±e/2). For one hole case, we study 
how the defect-defect correlation changes with t/J, which is the ratio between the hopping integral 
and cyclic exchange, equals to V/2t when 7>t Moreover, we argue that these fractional defects 
may propagate independently when the concentration of holes (or defects) is large enough. 

PACS numbers: 71.10.Fd, 71.27. +a 



I. INTRODUCTION 

During the past decade materials experiencing geomet- 
ric frustration has been a topic of much interest, both on 
experimental and theoretical sides [l| . The rich behavior 
in such systems is due to a large ground-state degeneracy, 
which renders them highly unstable to perturbations. A 
well known 3D frustrated structure is the pyrochlore lat- 
tice which is the backbone structure shared by many 
physically realistic materials. One may think of checker- 
board lattices as a 2D analog of pyrochlore lattices. This 
is the reason why a checkerboard lattice is often consid- 
ered in theoretical studies instead of the more realistic 
pyrochlore structure. 

Up to now most of work has been devoted to an un- 
derstanding of the magnetic properties of such frustrated 
lattices, while the charge degrees of freedom is waiting 
for more attention. Actually, the charge degrees of free- 
dom may exhibit fascinating physical properties in frus- 
trated systems such as pyrochlore and checkerboard lat- 
tices [IS]. Recently, one of the fascinating predictions 
on charge degrees of freedom has been proposed, say, the 
geometric frustration present on the pyrochlore lattices 
may give rise to fractional charges in two or three dimen- 
sions based on strong nearest neighbor repulsion close 
to half-integer filling. This proposal comes from the so 
called " tetrahedron rule" which is firstly stated explicitly 
by Anderson 0] to explain the metal-insulator transition 
in the spinel Fe304 [g], where the observed entropy re- 
duction is much less than expected from electrons with- 
out the strong short range correlation. Both pyrochlore 
and checkerboard lattices are made up of corner sharing 
units which are tetrahedron or crisscrossed plaquettes. 
The nearest neighbor repulsion will be minimized when 
each of the corner sharing units contains only two parti- 
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cles. If the kinetic energy can be neglected, the system 
will possess a large ground-state degeneracy, (4/3)' 3Ar / 4 ) 
for checkerboard lattices Taking (putting) one par- 
ticle from (into) the system, two tetrahedra (plaquettes) 
will emerge, each of them contains an extra particle (or 
hole). We will call this kind of tetrahedra (plaquettes) 
"defects". A demonstration of the configuration subject 
to "tetrahedron rule" and the formation of "defects" is 
shown in FigfT] Provided that the perturbation such as 
kinetic energy can be neglected, these two defects induced 
by one particle (hole) will propagate independently. If a 
particle carries charge e, one of the defects will carry 
charge e/2. Thus, such defects are "fractional". Note 
that the above argument is valid for both hardcore bosons 
and fermions as well as spinless and spin-S 1 system. On 
the other hand, current laser cooling and cold atom tech- 
nique make it possible to realize these systems on some 
artificial optical lattices 

However, the virtual processes induced by kinetic en- 
ergy will lift the high degeneracy and lead to differ- 
ent ground states. Because of high degeneracy, any 
small perturbation may change the low energy states vi- 
olently. To study this kind of quantum effects in spinless 
fcrmionic system on checkerboard lattices, Ref. [3j used 
exact diagonalization (ED) technique up to 32 sites. As 
a related problem to hardcore bosons, XXZ Heisenberg 
model in the Ising limit was studied by ED on small size 
lattices too Q . The ground state was identified as a non- 
magnetic state of resonating square plaquettes. Due to 
the limit of small size, there is lack of direct evidence to 
address the issue of confinement or dcconfincmcnt. 

In this paper, we will study spinless hardcore bosons by 
Green function Monte Carlo (GFMC) method which can 
give precise results on larger lattices (up to 28 x 28). The 
ground state energy is obtained with high accuracy and 
defect-defect correlation is present, which makes sense 
only on larger lattices. Also, we will inspect the issue 
of charge order. The outline of the paper follows. In 
section II, we introduce the extended Hubbard model and 
derive its effective Hamiltonian subjected to tetrahedron 
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rule in the strong repulsion limit. The relations between 
this model and XXZ Heisenberg model is discussed. In 
section III, a brief discussion of GFMC method is present, 
we also compare some results on small lattices with exact 
diagonalization. Section IV contains the main numerical 
results. The final section is devoted to conclusions. 



II. THE EFFECTIVE HAMILTONIAN 

We consider strong on-site repulsion U and nearest- 
neighbor repulsion V between spinless hardcore bosons 
on checkerboard lattices, where the intra-site interaction 
has been assumed to be a higher energy scale and ignored. 
Then the extended Hubbard Hamiltonian is of the form 

H = ~t + h.c.) + + U^2 n l C 1 ) 

(ij) (ij) > 

where bi(b\) denote annihilation (creation) operators at 
site i, rii = b\bi, and (ij) refers to a pair of nearest neigh- 
bor. We shall assume that U ^> V, t is large enough to 
ensure the non-double occupancy condition, only empty 
and singly occupied sites ar considered. Then we will fo- 
cus on the strong interaction regime V t > around 
half filling where the average occupation number per site 
is (rii) = 1/2, more precisely, half of the sites are occu- 
pied and the other half of the sites are empty. In that 
case the tetrahedron rule is imposed by the strong nearest 
neighbor repulsion, 

X> = 2, (2) 

where i belongs to a same crisscrossed plaquette. In this 
way, the inter-site repulsion will be minimized. Thus, 
in the limit V S> t, we obtain the effective low energy 
Hamiltonian in the subspace restricted by the condition 

Hj = -jY,b\b i2 b\b^ + h.c, (3) 

□ 

where J = 2t 2 /V, and □ denotes a four-site loop without 
crisscross, formed by sites i^i^ii- Away from half filling, 
we consider the hole doping only, due to the particle-hole 
symmetry at half-filling. Then the tetrahedron rule has 
to be modified as the following, 

and the effective Hamiltonian in this subspace is given 
by 

H eff = -tY J {b\b j + h.c.)+H J . (5) 

m 

It means that if we take one particle from the half filling 
system, there will emerge two defects where X^ieK Ui = ^ 



and Xieix] n i — 2 elsewhere. The hopping terms (propor- 
tional to t) of ([5]) will change the defects position but Hj 
will not. Although t/J = V/2t ^> 1/2, we can generally 
consider an effective model in which the ratio t/J ranges 
from zero to positive infinity. As pointed out in the for- 
mer references @, H, Q , if the virtual process at order J 
can be neglected, these two defects will propagate as in- 
dependent fractional objects. However, quantum effects 
such as the cyclic exchange of the order of J may or may 
not confine these fractional defects. 




FIG. 1: Left: A typical configuration subject to the tetrahe- 
dron rule at half filling. Right: Two defects induced by one 
hole may be considered as point like fractional excitations, 
whose motion is driven by hopping Hamiltonian. 

A hardcore boson model can be mapped onto a spin- 
1/2 model in general. The occupation number m cor- 
responds to the spin component Sf + |, and creation 
(annihilation) operators b\ (bi) correspond to S^(S^). 
The according spin- 1/2 model of (JTJ) is an XXZ Heisen- 
berg model. It is noted that this XXZ Hamiltonian is 
different from an usual XXZ model by the minus sign 
before spin flip terms, which can not be mapped to a 
positive one by an unitary transformation on a frustrated 
lattices. The related effective Hamiltonian in the Ising 
limit is a cyclic exchange which was studied by Shannon 
et al. [8j through ED on small lattices. The authors ar- 
gued that spinons (defects) in the above XXZ model will 
be confined due to the ground state correlation @. In 
this paper we confirm this argument through calculating 
the defect-defect correlation directly by GFMC method 
up to 24 x 24 lattices which are large enough to suppress 
the size effect. If only one hole is present the two induced 
defects will be confined together. However, as concentra- 
tion of holes and the ratio t/J — V/2t increase these 
defects may behave as independent point like objects. 



III. GROUND STATE WAVE FUNCTION AND 
GFMC APPROACH 

Note that both Hamiltonians © and ([5]) have only 
non-positive off-diagonal matrix elements in the Fock 
representation \nnni2 • • • } ■ It implies that the many- 
body boson wave function can be chosen to be non- 
negative everywhere in the ground state. This prop- 
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erty of the wave function will be crucial in applying the 
GFMC method to this system. 

GFMC is a general scheme for finding the lowest eigen- 
value of an operator. A trial eigenvector is subjected to 
a time evolution whose effect is to enhance those com- 
ponents of the solution with the lower eigenvalues of the 
operator. The ground state energy can be obtained as 
a mixed estimate. We choose the trial state \ipr) as 
an equal amplitude superposition in an invariant sub- 
space of Hamiltonian which contains the ground state. 
For a local operator O whose matrix elements satisfy 
(R\0\R') = S(R- R')O(R), such as density-density 
correlation function, which does not share eigenstates 
with the Hamiltonian, we use the " forward- walking" 
technique, well known in many-body theory [l(| [ll|, [lj, 
[T3I fl4| . to measure the expectation values. In this way, 
at the nth step of iteration, O (R^) is evaluated for each 
configuration R} n in the ensemble {R n }- For the remain- 
ing to steps of the random walk, a record is kept of the 
configuration of {-R^} from which each subsequent en- 
semble member has evolved. At the end of the m + nth 
step, (O) is evaluated by forming the weighted summa- 
tion [n[ 

(O) = — ^ J -, 

Ej V>T [Rin+n) 

where the notation i (J) indicates that i is the progenitor 
of j. In this paper, we will calculate density-density cor- 
relation and defect-defect correlation, both of them are 
local operators and this method can be applied. 



A. Comparison with exact diagonalization 

GFMC is a very accurate method to obtain the ground 
state and low-lying excited state property of many-body 
interacting systems. In the past it has been applied 
successfully to the ground state properties of helium, 
interacting electron gas, small molecules, Heisenberg 
model on 2D square lattices, and lattice gauge theories 
0, HO, El GJ, m- To check the precision of the GFMC 
method in our system, we will compare some results from 
the ED with those from GFMC at small lattices with pe- 
riodic boundary condition (PBC). Hereafter we will use 
PBC in this paper. 

The ground state energy at half filling and at one hole 
doped lattice are given in Tables U and [H] respectively. 
The digits in the bracket are statistics error bars. 



TABLE I: The ground state energy E/J at half filling. 



Size 


4x4 


4x6 


6x6 


ED 


-4.47214 


-6.46995 


-9.47393 


GFMC 


-4.4717(12) 


-6.4695(19) 


-9.4739(22) 



TABLE II: The ground state energy E/J of one hole doped 
4x4 lattice. 



t/J 


ED 


GFMC 


0.02 


-3.52850 


-3.5290(9) 


0.2 


-4.91215 


-4.9119(16) 


2.0 


-22.92923 


-22.9310(16) 



Density-density correlation D(i,j) = D(i — j) = 
(niTij) — j at half filling on 6 x 6 lattice have been exam- 
ined too. The results from ED and GFMC are present 
in Table Mil and IIVI respectively. Each of the two ta- 
bles contains a matrix D(i x ,i y ) whose column and row 
indices i x , i y = 0, 1, 2, 3 correspond to the displacements 
along two directions respectively. 

TABLE III: Density-density correlatoin D(i x ,i y ) on 6 x 6 lat- 
tice calculated by ED. 



0.25 


-0.13170 


0.01933 


-0.02527 


-0.13170 


0.07083 


-0.01590 


0.02185 


0.01933 


-0.01590 


0.01584 


-0.01920 


-0.02527 


0.02185 


-0.01920 


0.01998 



TABLE IV: Density-density correlatoin D(i x ,i y ) on 6 x 6 lat- 
tice calculated by GFMC. 



0.25 


-0.1319(7) 


0.0196(12) 


-0.0254(11) 


-0.1319(7) 


0.0711(9) 


-0.0161(10) 


0.0220(8) 


0.0196(12) 


-0.0161(10) 


0.0160(10) 


-0.0193(7) 


-0.0254(11) 


0.0220(8) 


-0.0193(7) 


0.0200(7) 



The staggered pseudospin magnetization (|7|) on a half 
filled 6x6 lattice given by GFMC is 0.1038(8), while the 
result from ED is 0.10417. 

From the above, one sees that GFMC algorithm is an 
effective method to deal with spinless hardcore bosons on 
checkerboard lattices. 



IV. RESULTS 

Simulations were carried out for L x L lattices up to 
L = 28 (L = 24 for correlations). Time steps from 0.003 
to 0.05 were used, depending on lattice size and the ratio 
t/J, here we set J = 1. In practice, we use about 100L 2 
generations to reach the ground state distribution, then 
iterate about 1000L 2 generations to measure the physi- 
cal quantities. To avoid self-correlation and improve the 
efficiency, we make an expectation value measurement 
only after every L 2 iterations. It is not advantageous 
to perform measurements at still larger intervals, since 
the measured generations are already nearly statistically 
independent. For ergodicity we should control the pop- 
ulation of random walkers large enough. Otherwise they 



4 



will be trapped in a higher energy state instead of the 
ground state. According to our experience, keeping L 3 
random walkers in each generation is enough to ensure 
ergodicity on L x L lattice. 



A. Half filling: ground state energy and staggered 
pseudospin magnetization 

Firstly, we calculate the ground state energy at half 
filling which may serve as a standard to compare with 
other analytical or numerical study in this system. Fig. 
[5]shows ground state energies per site up to 28 x 28 lattice. 
We extrapolated it to the thermodynamic limit through 
the following formula, 

E{L) _ E x E 2 (I 

-r- Ea + - + L^ + {-^ 

where E = -0.2591(4), E x = 0.008(9), E 2 
and L is the linear size of lattices. 
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FIG. 2: Ground state energy per site at half filling and ex- 
trapolation to thermodynamic limit. Values are shown for 
L x L lattices, L = 8, 10, • • • , 28. 



because if we map a hardcore boson model to a spin-1/2 
model this quantity is nothing but the staggered magneti- 
zation. If the system is charge ordered, m will not vanish 
in the thermodynamic limit. By "forward walking" tech- 
nique, we calculate the pseudospin magnetization up to 
24 x 24 lattice and extrapolated it to the thermodynamic 
limit by the following formula, 



rn i 



m (L) = m + — + -j2 + ^ 



(8) 



with m = 0.000(4), mi = 0.55(9) and m 2 = 0.6(5). 
As shown in Fig[31 it results in a vanishing charge order 
at (7r, 7r) in the thermodynamic limit. This result agrees 
with the claim of non-magnetic phase in Ref. Q . Similar 
conclusion was found for spinless fermion Q. 
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FIG. 3: Staggered pseudospin magnetization, Eq. (O, on 
finite size lattices at half filling and extrapolation to ther- 
modynamic limit. Values are shown for L x L lattices, 
L = 8, 10, • ■ ■ , 24. Zero result for mo means charge order 
is absent here. 



Another issue that we shall address is whether there ex- 
ists charge density order in the thermodynamic limit. Al- 
though density-density correlation D (i,j) exhibits (7r, tt) 
charge-modulation patterns on small lattices with PBC, 
as it behaves in fermionic system |3( , it may vanish in the 
thermodynamic limit. To answer this question, we intro- 
duce the staggered pseudospin magnetization m, which 
is defined as 




where N is the number of sites on checkerboard lattices. 
The staggered pseudospin magnetization has its name 



B. One hole doping: ground state energy and 
defect-defect correlation 

One of the central questions is whether or not the de- 
fects will be confined by the Hamiltonian (0 . To answer 
this question, we calculate the ground state energy at 
one hole doping and t = at first. When t = 0, the two 
defects can not move away by hopping terms, so that 
we should calculate the ground state energy with fixed 
defects. Fig. [4] shows the numerical results for L x L 
lattice, L = 16, 20, 24. It turns out that the ground state 
energy will increase linearly as the distance between two 
defects increases. The distance between two defects is 
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defined as the distance from one plaquette center to an- 
other in units of the lattice constant a, say, it is y/2 for 
two nearest neighbor plaquettes. 




4 6 8 10 12 14 16 
Distance between two defects (r) 



20 



FIG. 4: (Color online) Cyclic exchange energy Hj will in- 
crease linearly when the distance between two defects in- 
creases in the case of t = 0. Here Eq is the ground state 
energy at haif fifling. 
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FIG. 5: (Color online) The defect-defect correlation on 20 x 20 
lattice. As the ratio t/J increasing, it will become more and 
more uniform. Inset: The average distance R between the 
two defects increases monotonically with t/J. The saturation 
is due to the size effect, for a homogeneous distribution, the 
maximum R wiil be 0.3826L. Sofid dots are piotted for 20 x 20 
fattice and circies for 16 x 16 lattice. 



Even though we turn on the hopping term, the upper 
bound of the gained kinetic energy is less than 2zt, z = 
6 for checkerboard lattices is the coordination number. 
Since the confinement potential increases linearly with 
the distance between the two defects, the two defects 
will be confined with an average distance R. 

Now we turn to t > case and study how the defect- 
defect correlation will change as the ratio t/J varies. The 
defect-defect correlation C (p, q) = C (p — q) can be de- 
fined as: 

C(p,g) = /(2-53 ni )(2-X;n J -A, (9) 

\ iep jeq I 

where p, q denote the crisscrossed plaquettes. The nu- 
merical result is present in FigO One notices that when 
t/J is small, the two defects will be linearly bounded. As 
t/J increasing, the defect-defect correlation will be more 
and more uniform in a finite size lattice with periodic 
boundary condition. Hence we expect that the average 
distance R(t/J) between the two defects in the infinite 
lattice will increase monotonically with t/J. However, as 
shown in the inset of FigJSJ it will reach the saturation 
on a finite lattice due to the size effect, the maximum 
will be (V2 + ln(l + V2))L/6 = 0.3826L on L x L lattice 
for a homogeneous distribution. 

From the above we can conclude that the cyclic ex- 
change Hj provides a linear confinement potential be- 
tween two defects but the hopping term tends to increase 



this average distance. If the concentration of defects 
S = 2Nhoie/N is large enough as to satisfy 

when the average area occupied by one defect, 
No 2 / (2Nhoie), is smaller than the confinement area, ttR 2 , 
the defects will have a homogeneous distribution on the 
infinity lattice instead of to be confined together in cou- 
ples. It implies that it is possible to treat these defects 
as independent point like excitations. 

V. CONCLUSION 

In summary, we apply the GFMC algorithm to a spin- 
less hardcore boson model with strong nearest neighbor 
repulsion on checkerboard lattices near half filling. To 
avoid any variational bias from the trial wave function, a 
"forward walking" technique has been used to compute 
density-density correlation and defect-defect correlations. 
At half filling, the ground state energy is obtained and 
extrapolated to infinity size. It turns out that the stag- 
gered pseudospin magnetization vanishes in the thermo- 
dynamic limit, thus charge order at (n, n) is absent in this 
system. Away from half filling, two defects induced by 
each hole (particle) may carry fractional charge (±e/2). 
In the case of one hole doping, we study how the defect- 
defect correlation changes with the parameter t/J, which 
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equals to V/2t when V 3> t. The cyclic exchange is found 
to provide a linear confinement potential between two de- 
fects, while the hopping term as kinetic energy is tending 
to separate them away. Moreover, we argue that these 
defects may propagate independently when the concen- 
tration of holes (or defects) is large enough. 

Finally, confinement (deconfmement) is a subtle issue, 
although we have presented some evidences for the pos- 
sibility of the existence of defects which may carry frac- 
tional charge (±e/2) in this system, a lots of work remain 
to be done. This issue should be treated by other analyti- 



cal and numerical methods. Especially, the effective field 
theories for the present system are expected to describe 
the ground state and low lying excitations well, thus a 
confident conclusion for confmement/deconfincmcnt will 
be achieved. 
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